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ABSTRACT 


The  purpose  of  this  study  is  to  investigate  the  asymptotic  behaviour 
of  certain  (conditional)  rank  tests  for  comparative  experiments  proposed 
by  Hodges  and  Lehmann  (1962)  [8],  and  their  extensions  and  variations. 

The  results  obtained  in  this  thesis  are  extensions  of  the  results  of  Mehra 
and  Sarangi  (1965)  [ lU ]  (where  the  case  of  complete  designs  was  discussed) 
to  the  case  of  balanced  and  partially  balanced  incomplete  block  designs. 
The  methods  employed  are  based  on  certain  properties  of  the  U-statistics 
derived  by  Hoeffding  (19^8)  [9], 

In  Chapter  I,  a  brief  resume  of  the  relevant  earlier  work  has  been 
given  and  the  test  proposed  by  Hodges  and  Lehmann  described.  In 
Chapter  II,  a  simplified  large  sample  version  of  this  test  is  discussed, 
which  reduces  its  application  to  reading  the  critical  value  from  the 
chi-square  tables. 

In  Chapter  III,  the  asymptotic  distribution  of  the  test  statistic, 
under  a  sequence  of  alternatives  approaching  the  null  hypothesis  has 
been  found  to  follow  a  non-central  chi-square  law.  This  helps  us  to 
study  the  power  properties  and  the  Pitman  efficiency  of  the  test  statis¬ 
tics  in  Chapter  IV.  Explicit  numerical  evaluations  of  the  expressions 
for  asymptotic  efficiency,  which  has  been  limited  to  the  normal  form  of 
the  parent  distribution,  indicates  a  rather  high  degree  of  efficiency 
of  the  test. 

Chapters  IV  and  V  contain  several  remarks  that  follow  from  the 
results  of  our  investigation,  including  a  useful  lower  bound  for  the 
asymptotic  efficiency  of  the  investigated  test  relative  to  the  classical 
F-test.  Further  investigations  that  will  throw  more  light  on  our 
present  findings  have  been  suggested  at  suitable  places. 
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CHAPTER  I 


THE  PROBLEM  AND  THE  DESCRIPTION  OF  THE  TEST  PROCEDURES 

1.1.  Introduction 

The  emphasis  on  many  of  the  classical  test  procedures  has  consid¬ 
erably  diminished  in  the  last  decade,  in  relation  to  a  certain  class  of 
statistical  hypotheses,  with  attempts  to  find  distribution- free  tests. 
Almost  all  classical  tests  are  based  on  the  underlying  model  of  "random 
sampling  from  a  normal  population."  These  test  procedures  are  satisfactory 
for  most  practical  purposes  if  the  sample  is  "large"  and  the  sampling 
distribution  of  the  test  statistic  can  be  regarded  as  approximately  normal. 
But  it  has  now  been  realised  that  some  of  these  tests  are  vulnerable  to 
gross  errors  if  conditions  of  normality  are  violated.  Their  levels  of 
significance  are  also  not  exact,  unless  the  parent  population  is  normal. 

Non-parametric  methods,  which  overcame  the  above  defects  were 
introduced  quite  early  in  statistical  literature.  But  no  systematic 
study  of  this  branch  was  made  till  the  19^0 ' s .  In  the  last  twenty  years, 
this  area  has  been  subjected  to  intensive  study  and  many  new  non-parametric 
tests  have  been  proposed  and  their  proporties  investigated.  These  tests 
have  invariably  been  shown  to  be  preferable  to  the  corresponding  classical 
tests . 

The  scope  of  our  present  study  is  confined  to  the  investigation  of 
the  asymptotic  behaviour  of  certain  type  of  rank  tests  for  comparative 
experiments,  proposed  by  Hodges  and  Lehmann'  [8]  and  their  extensions  and 
variations.  Subsequently,  Lehmann  developed  a  new  class  of  non- 


2 


parametric  tests,  whose  theory  runs  parallel  to  that  of  normal  theory 
testing,  estimation  and  multiple  comparison  procedures.  However,  the 
efficiency  of  their  earlier  conditional  rank  test  has  remained  uninves¬ 
tigated  o 

Suppose  we  have  k  objects  or  "treatments"  (e.g.  varieties  of 
wheat,  agricultural  fertilizers,  products  from  similar  manufacturing 
units  etc.)  to  be  compared.  In  many  practical  situations  our  decisions 
have  often  to  be  based  on  experiments  conducted  on  certain  "experimental 
units."  However,  in  general,  different  experimental  units  may  react 
differently  towards  the  same  treatment  and  thus  any  method  of  precise 
comparison  should  divide  the  experimental  units  into  relatively  homo¬ 
genous  groups  called  "blocks."  Then  the  experimental  layout  can  be 
represented  as 
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where  there  are  m.  units  in  the  i • 0  block  corresponding  to  the  j** 
treatment  (or  in  the  (i, j )•£  cell).  In  what  follows,  asymptotic  means 

n 


"as  n  tends  to  infinity. 


( i )  ONE-WAY  LAYOUT 
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If  the  reactions  of  the  experimental  units  in  different  blocks 
towards  the  same  treatment  are  not  different,  then  our  analysis  reduces 
to  that  of  the  one-way  layout  or  the  k-sample  problem.  This  problem, 
in  the  special  case  when  k  =  2,  was  considered  by  Wilcoxon  [20]  in  1945. 
In  1948,  Pitman  [l6],  who  developed  the  concept  of  efficiency,  found  that 
the  efficiency  of  this  test  relative  to  the  corresponding  normal  theory 
test  (t-test)  is  (3/n  )  if  the  parent  population  is  normal.  In  1952, 
Kruskal-Wallis  [10]  generalised  Wilcoxon* s  test  to  any  positive  integral 
k;  and  the  corresponding  asymptotic  efficiency  was  also  shown  to  be  (3/n  ) 
(Andrews  [  1]  ,  195*0. 

(ii)  TWO-WAY  LAYOUT:  SEPARATE  RANKING 


However,  if  the  reactions  of  the  experimental  units  differ 

from  block  to  block,  the  design  reduces  to  a  two-way  layout „  The 

earliest  non--parametric  test  for  this  was  suggested  by  Wilcoxon  in 

1946  for  the  paired  comparison  case  when  k  =  2  and  all  m  =  1*  Pitman 

(1948)  found  the  asymptotic  efficiency  of  this  test  relative  to  the 

t-test  to  be  3 /II ^  if  the  parent  population  is  normal.  Friedman  [  5 ]  »  in 

1937 ,  proposed  a  test  based  on  independent  rankings  within  each  block  for 

the  case  when  all  m. .=  1  ice.  there  is  one  observation  per  cell.  Durbin 

i  j 

[4]  in  1951 s  extended  the  test  to  cover  balanced  incomplete  block  designs 
when  m.  .=  o  or  1.  The  asymptotic  efficiency  of  these  tests  was  computed 
by  van-Elteren  and  Noether  (1959,  [15])  to  be  {3m/n(m+l) } ,  where  m  is  the 
block  size.  In  1953,  Benard  and  van-Elteren  [2]  discussed  these  tests  in 


great  detail  and  proposed  tests  for  a  wide  class  of  incomplete  block 
designs.  Hodges-Lehmann  [7]  showed  that  the  asymptotic  efficiency  of 
Wilcoxon  test  relative  to  the  t-test  never  falls  below  .864,  although  its 
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value  can  be  arbitrarily  high  and  even  infinite.  This  very  desirable 
feature  of  attaining  a  high  efficiency  independently  of  what  the  parent 
population  is,  has  been  called  "robustness"  in  statistical  literature. 


TWO-WAY  LAYOUT:  JOINT  RANKING 


It  is  noted  that  for  k  =  2  and  all  m.  =  1,  the  paired  com- 
parison  test  attains  an  asymptotic  efficiency  of  3/n  relative  to  the  t- 
test  against  2/n  of  Friedman's  test.  The  higher  effiency  of  the  former 
seems  to  be  due  to  the  fact  that  it  takes  account  of  certain  interblock 
comparisons  which  are  entirely  ignored  by  the  latter.  A  first  attempt 
at  the  investigation  of  a  test  statistic  based  on " interblock  comparisons 
was  done  by  Mehra  [ 13 ]  for  m  =  2  and  general  k;  and  the  corresponding 
asymptotic  efficiencies  agreed  with  that  of  the  paired-comparison  test 
relative  to  the  competing  test  procedures.  Mehra  also  extended  the 
results  to  certain  paired  comparison  experimental  designs.  Hodges  and 
Lehmann  [8],  in  1962,  also  proposed  a  conditional  test  based  on  a 
combined  ranking  of  all  the  observations  which  thus  takes  into  account 
some  interblock  comparisons  of  cell  observations.  In  this  chapter,  we 
discuss  an  asymptotic  version  of  the  proposed  test.  The  asymptotic 
efficiency  of  the  case  of  complete  experimental  designs  was  discussed 
by  Mehra  and  Sarangi  [lU] .  Here  the  results  of  Hodges-Lehmann  [8]  and 
Mehra-Sarangi  [ lU ]  are  extended  to  other  general  designs  like  the  balanced 
and  partially  balanced  incomplete  block  designs  etc. 

1. 2,  Notations 


Let  X, 


ij 


(  i  1 1  2 , 


,n;j  =  1,2,  • . ,k;l  =  1,2, 


»m. 


ij 


)be  the  ^ 


th . 


observation  in  the  (i,j)th.  cell.  We  assume  these  m_  independent  random 
variables  (r0v0)  to  be  identically  distributed  with  a  common  continuous 


-  5  - 

cumulative  distribution  function  (c.d.f.)  F..(x)  such  that 

^  J 


(1.2,1) 


to  '•  F  (x)  =  F  (x  +  C  .) 

ij  J  1 


where  S^s  are  unknown  constants  called  "block  effects."  The  null 
hypothesis  of  no  treatment  differences  is  equivalent  to 


H 


Fi  =  F2  = 


=  F, 


o  i  **  k 

We  note  that  model  (1.2,1)  is  the  assumption  of  additivity  of  block 

effects  in  the  classical  linear  analysis  of  variance  modelo 

1 o 3 «  Method  of  Alignment 

The  conditional  test  proposed  by  Hodges-Lehmann  [8]  is  based  on 
a  combined  ranking  of  all  observations  in  different  blocks  after  they 
have  been  made  comparable  by  removing  the  block  effects.  The  removal  of 
the  block  effects  is  referred  to  as  "alignment."  One  method  of  alignment 
is  to  subtract  from  each  observation  in  a  block,  say  the  i"^ ,  a  function 
y  of  the  observations  in  that  block,  satisfying  the  condition. 


{U3‘1)  p(xill  +  a-  xii?  +  a-  •••>  xilm.  + 

'  ll 


a,  .  .  ,x..  +  a ,  .  .  .  ,x „ +  a) 

lkl  lkm. 


y (*il 1  ,  xi 1 2 » *  0  *  >xil 


m 


»xiki » •  •  •  ‘X-i 


il 


ikm , 


ik 

)  4  a 


It  may  be  observed  that  the  aligned  observations  in  each  block  should 
preferably  have  a  joint  symmetric  distribution  in  its  arguments  under  the 
hypothesis 0  This  condition  is  not  necessary  e.g.  our  method  of  alignment 
may  be  such  that  the  joint  distribution  is  symmetric  in  disjoint  subsets 
of  its  arguments.  But  we  confine  our  attention  only  to  the  former  case, 
which  is  satisfied  if,  in  particular,  y  is  a  symmetric  function  of  its 
arguments t  This  is  shown  by  the  following  lemma. 


LEMMA.  1,1,  Let  Xi,X2 


3  o  o  o  3 


X  be  independent  identically  distributed  r.v.’s 


and  let  y(X1  X?  „ „ „ c ,X_)  be  a  symmetric  function  of  these  variables.  Then, 

*  j  c 

the  r.v.'s  Zi  Zo  0 . 0 „  ,Z  where  Z.  =  X.  -  y(Xi  X2  . . . . ,X  )  have  a  symmetric 
-  »  ,  c  — ■ —  j  j  ’  5  c  - 
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joint  distribut i on . 

Proofs  This  lemma  was  proved  in  [lU]. 
F(t2}ti}Ot0C}tc)  =  Pr[Zj  <_t2^Z2 

=  Pr [Xx  ~  y (X !  X2  co, 

>  * 


We  also  give  the  proof  below. 


0  0  f  O  ,  Z  <S  t  J 
c  —  c 

5XC  )  <_  t2  X2  ~  y(Xi  X2> 


,X) 


<  ti  Z3  <  t3  .  .  .  •  ,  Z  <  "C  ] 

>—>  1  >  J -  0  5  5  C  _  c 


=  Pr[X2 


t(x25xl5 


,x  )  <  to,X 


c  — • 


-  y(X2,X1,...,Xc) 


=  Pr[Z2 


t2  Z,  c  ti  ...,Z  <  t  ] 

—  5  C  — •  C 


=  F(t! }t2j . . . ,tc) 
by  the  symmetry  of  u(Xi  sX2  ,  .  .  .  ,X  ). 

It  is  seen  that  statistics  like  the  arithmetic  mean,  median  etc.  satisfy 

the  above  conditions, 

lebo  Description  of  the  test 

A  fairly  complete  discussion  of  the  test  and  its  asymptotic 

versions  are  given  in  [8]  for  the  special  case  of  k  =  2.  Our  discussion 

concerns  the  asymptotic  version  of  the  test  for  general  k.  Let 

Z  (i  =  l,2,...,n;j  =  1,2, ...k;  j,  =  1,2,. .  .  ,m  )  denote  the  aligned 
iJi  ^ 

observations  and  r. .  denote  the  rank  of  Z . .  in  a  combined  ranking 

ijl  iJl 

of  all  the  N  =  E£m„  .  aligned  observations.  Following  [8],  each 

ij  ij 

conditional  event  when  the  set  of  ranks  in  each  block  is  given  will  be 
referred  to  as  a  "configuration.  M  Given  a  configuration,  under  the 
null  hypothesis  of  no  treatment  differences,  the  ranks  within  each 
block  can  be  assigned  to  the  'treatments  at  random,  since  the  vector  of 
aligned  observations  within  each  block  has  a  symmetric  joint  distribution. 
The  following  discussion  assumes  (i)  that  the  method  of  alignment  is  the 
same  in  each  block  and  (ii)  that  conditions  (A)  and  (B)  of  [8],  given 


below  are  satisfied. 
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(A)  After  alignment,  each  block  contains  at  least  one  positive  and 
one  negative  observation. 

(B)  Complete  randomization  is  used  within  each  block. 

These  conditions  are  very  general  and  satisfied  for  most  methods  of 

alignment  like  those  on  the  mean,  median  or  the  Winsorized  and  censored 

means.  Let  N.  =  £m.  .  R.  =  ZZr .  . ,  =  sum  of  the  ranks  under  the 
i  J  iJj  J  1 1  iJ  l 

treatment.  Following  [lU],  let  E(*)»  var(*)  and  cov( • , • )  denote  the 

conditional  expectation,  variance  and  covariance  respectively,  under 

H  ,  given  a  configuration.  Then  by  direct  computation,  we  get 

n 


{l.k.D 


E(R.)  =  Z  m . .  r. 

J  i=1  U  1 


2  n 

a.  =  a  =  Var(R.)  =  Z  (N -  m .  . )  m ..  {T?  /  (N.-  1)} 

j  jj  J  i=l  1  ij  i  i  ’ 

n 

a  t  =  Cov  (R .  ,R .  i )  =  -  Z  n.  m.  .  ,  {T.2/(N.-  l)} 

JJ  J*  J  i=1  iJ  iJ  i  i 

where  T.2  =  -££N;( r .  -  r.  )2  and  r.  =  {ZZ  r.  /N .  ) .  We  now  prove: 

i  ij  t  i  l  .  iJl  i 


Jl 


Theorem. 

If  ( c ' )  every  pair  of  treatments  occurs  together  at  least  in  one 
block,  then  the  variance-covariance  matrix  has  rank  ( k- 1 ) . 

Proof:  The  theorem  is  a  direct  consequence  of  the  following  result  due 

to  Taussky  [l8]. 

Theorem  (A).  If  (a..*)  be  a  kxk  matrix  of  complex  elements  such  that 

J  J 


a 


JJ 


k 

>  ^  I  0  i  1 1 
J  T#J  =  1 


(j  =  1  ,2,  .... ,k) 


then 


(ojj,)|?* 0 


In  our  case ,  o . . > 0,  a . . »  <  0  and  a .  .  =  £  a • • f  •  Thus  the  matrix 

JJ  JJ  Jj  j  =  i  JJ 

detained  from  (a,  . »  )  by  deleting  the  row  and  column  satisfies 

J  J 
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the  conditions  of  the  theorem,,  Hence  this  matrix  and  thus  the 
variance-covariance  matrix  (singular)  has  rank  (k  -  l). 

The  proposed  conditional  test  statistic  of  [8]  can  now  he  obtained  as 
follows 6  Consider  the  matrix  obtained  from 


all 

a,  o 

i  Z  c*oooooe 

alk 

■f  R 1 — E  ( R  i ) } 

a12 

p  oooooooo 

a2k 

(r2-e(r2)} 

Oik 

a2k  ©ooooooo 

°kk 

{vi(v} 

(Ri-ECRi)}  {R2-E(R2)}  o,.  (Rk-E(Rk)}  o 

by  deleting  the  row  and  the  column  (jfk+l)  and  compute  its 
determinant  |AU|;  consider  the  above  matrix  without  the  last  row  and 
last  column  and  compute  the  determinant  |a|  of  the  matrix  obtained  by 
deleting  the  same  row  and  column  from  it.  The  test  statistic  is 
(I.U.2)  Wn=|Au|/|A| 

The  test  consists  in  rejecting  Hq  at  a  level  of  significane  a  ,  if 

exceeds  a  predetermined  number  a  such  that  the  conditional  probability 

PH  [Wn  >  =  a»  The  application  of  this  test  requires  the  null 

distribution  of  W  which  can  be  done  by  random  allocation  of  the  ranks 

n 

within  each  block  to  the  treatments  present  there.  A  computing  machine 

may  be  used  unless  the  design  is  small  for  manual  operations. 

The  statistic  W  takes  a  more  convenient  form  when  m . .  does  not 

n  lj 

vary  from  block  to  block,  viz. 

(2.4.3)  m .  .  =  m2  •  „  o  o  o  o  —  m  .  —  m .  1  °  j  (j  —  •••?  k) 

lj  J  J 

We  give  here  details  of  this  case,  which  was  dealt  with  in  [l4].  Under 

(2. 4.3)  9  the  expressions  (2.4.2)  simplify  to 


-  9  - 


E(R  )  =  m , n ( nN ’  +  l)/2 
J  J 

a  =  Var (R  )  =  { (N*  -  m  )m  / (N '  -  l)}  (  hj  ) 
J  J  J  J  J  i =1 

Ojj ,  =  CovfRj.Ej,)  =  -  {mJmJ,/(N'  -  l)}  (IT?) 


k 

where  N'  =  Em,  Here,  the  inverse  of  the  variance-covariance  matrix, 

j=l  J 

after  deleting  the  last  row  and  last  column  is  easily  verified  to  be 

A  1  =  ( e  .  .  ,  ) ,  where , 

J  J 


N 


1 


n  n ' (  t  t? ) 


/ 

l  m. 


k 


1_ 

m 


i=l 


ijj,  =  +  (N*  -  1)/N'mk( 


n 

Z 

i=l 


T2) 

l 


(j  i  y  =  1,2,.. . ,(k-l)) 


It  is  easily  verified,  using  the  relation  Z  {R 


m  n(nl' 
J 


+  l)/2}  =  0 


that , 


(l.k'k) 


W  = 
n 


(N'  - 


n 

1)/N'(  Z  T2) 

1=1  "J 


k 

Z 


_1 

m 


(R  .  - 


J=1  J 


m  c  n  ( nN  * 
J 


+  l)/2}: 


For  the  simpler  case  of  one  observation  per  cell,  which  has  drawn  special 

attention  in  statistical  literature,  we  get 

'  k  “ 

Z  {R.  -  n(nk  +  l)/2}2 

J=1  J 


W  = 
n 


(k  -  l)/k( 


n 

Z 

i=l 


Tj) 


CHAPTER  II 


ASYMPTOTIC  VERSION  OF  THE  TEST  UNDER  THE  HYPOTHESIS 


2 c 1 o  Asymptotic  Null  Distribution 

It  has  been  pointed  out  in  the  last  chapter  that  the 

computation  of  the  null  distribution  of  W  is  rather  cumbersome  and 

n 

the  magnitude  of  labour  increases  rapidly  with  increase  in  sample 
size,  However,  when  the  sample  size  or  the  number  of  blocks  is 
large,  the  critical  value  ^  can  be  approximated  by  using 

V2 

x  -  distribution  tables,  as  shown  by  the  following  theorem.  It 
is  important  to  note  that  in  the  following  theorem,  the  convergence 


of  to  the  X  _  variable  is  uniform  in  the  configuration 


Theorem  2,1,  Suppose  the  hypothesis  Hq  is  true,  the  method 
of  alignment  be  such  that  (A)  and  (B)  hold  and  also  the  conditions 

(C)  as  n  tends  to  infinity,  {N„.,(n)/n}  is  bounded,  awav  from  zero 

J  J 

for  each  pair  ( j  ,  j  '  ) ,  where  N . . , (n)  denotes  the  number  of  blocks 

J  J 

in  which  both  the  treatments  T.  and  T „ „  occur  together  and 

‘ - — — — - ~ — - -  J  -  J  ’  -  - 

(D)  N.  <  n  for  all  i.  Then  (i)  the  conditional  test  statistic  W  , 

l  - - —  — — -  - - — — - - — - - - - -  n 

given  a  configuration,  converges  in  distribution,  as  n  tends  to 

infinity,  to  the  X2  -  variable  with  ( k— 1 )  degrees  of  freedom:  and 

( i i )  the  convergence  is  uniform  in  the  configuration. 

This  theorem  was  proved  by  Mehra  and  Sarangi  [l4]  under 

the  restriction  m_  >_  1  i,e,  for  complete  experimental  designs. 

These  results  are  extended  below  to  the  case  where  some  m, ,'s  are 

ij 


zero  i,ec  incomplete  block  designs. 


11 


k-1 


m.  . 

Let  R.  =  £  r. 
ij  1=1  iJ* 1 


Consider  an  arbitrary  linear  conbination 


U.  =  Z  c .R , .  of  R. J  =  1,2, . . . , (k  -  l) .  Let  b?  and  3 .  denote  the 
1  J  1J  1J  11 

j  O  n  9 

conditional  variance  and  third  absolute  moment  of  U.  and  let  Sz  =  £  bf. 


n 


i=l 


The  proof  of  the  theorem  follows  from  the  following  Lemmas . 

Lemma  2.1.  Under  the  assumptions  (A),  (B),  (C)  and  (D),  (i)  the 


n , 


conditional  (given  a  configuration)  distribution  of  E{U„  -  E(U.)}/S 

i=l  1  in 

converges,  as  n  -*«>,  to  the  standard  normal  distribution;  and  (ii)  the 

convergence  is  uniform  in  the  configuration. 


Proof o  In  the  sequel,  let  S.  denote  the  set  of  all  blocks  containing 
m.  J 

the  j"!^  treatment  and  S  ,  denote  the  set  of  all  blocks  containing  both 

<3  J 

the  o  and  treatments.  The  following  possibilities  can  occur. 

(i)  Suppose  c.  #  0  and  c.t  =0  for  all  j*  ^  j.  Then,  there  exists  a 

j  j  1 

constant  0  <  a^  <  °°,  independent  of  the  configuration  such  that 


n 


n 


(2.1.1)  £  3.  <  |  c  .  |  3  J  E  |  R  .  .  -  E  R.  ,  |  3  <  a^1 

L  i  -  1  J  1  L  1  ij  ij  1  -  1 

i=l  i=l 

Further,  from  (l.i+.l), 


n 


S2  =c2  I  (N.  -  m.  .  )m.  .{T?  /(N.  -  l)}  >  c  2 


V 


L  T_.  2/(n~l) 

isS «, 

J 


n 


i  L 
I 

i=l 


ij  1J  1 


J  ,  i 


By  assumption  (C),  there  exists  0  <  f  <_  1,  independent  of  j  such  that 

for  every  n  >  n^  (to  be  chosen  large  enough),  there  are  more  than  nf 

elements  in  S  .  If  the  original  aligned  observations  consist  of  s 
J 

negative  and  (N  -  s)  non-negative  numbers;  every  positive  observation 

must  have  a  rank  greater  than  s  and  every  negative  observation  a  rank 

not  greater  than  s.  In  any  case,  by  (A),  for  every  i  there  exists 

a  r . ,  such  that , 
ijl 
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r .  ...  -  r .  i  > 
ljl  i 


r.  -  s 
ijl 


Thus , 


V 

L 

ieS  „ 
J 


=  l  l  l  (rij1 

ieSj  J  1 


Th )  2/Ni  >_ 


L  ljl 


-  s 


)2/ 


n 


IeS. 

cJ 


>  2 


lz+  2Z+ 


+  [(nf  -  1 ) /2 ] 2 


/ n  >  a£n3 


where  0  <  a2  <  00  is  independent  of  the  configuration  and  [p]  denotes 
the  greatest  integer  contained  in  p;  since  under  the  continuity  of  the 
c.d.f.'s  there  will  he  no  ties  for  non-zero  aligned  observations  with 
probability  one.  Thus,  there  exists  0  <  a2  <  00 ,  independent  of  the 
configuration  such  that, 


(2.1.2) 


S2  >  a2n3 
n 


(ii)  Suppose  cx  =  c2  =  . . .  =  ?  0 .  Then 

n 

Z  U  =  c  [  N(N  +  l)/2  -  R  ] 

i=l  1  1  k 


since  the  sum  of  all  the  ranks  is  N(N  +  l)/2.  Equations  (2.1.1)  and 
(2 .1., 2)  can  now  be  proved  as  above 0 

(iii)  Suppose  there  exist  at  least  two  c’s,  say  c.  and  c  which  are 

j  y 

not  equal.  Then  (20lol)  follows  as  before  and  (2.1.2)  is  proved  below. 

n  n  r  k-1 

V 

L 


S2  =  y  b2  =  1 
n  L  i  L 


i=l 


l  c2  (N_.  -  m_j  .  )m_j  .  -  2  ll  c.c.,m.  .m_ 


i=l  L  j=l 
“  r  k-1 


L 

i=l 


ij  ij 


k-1 


i 


J  J  ij  ij 


x[t2/  (N\  -  1)] 
k-11  1  k-1 


1  ci  ^  z  Zc.m  )/(  z  m 

j=i  J  J  j=i  1J  j=ij  1J  j=i 


k-1 

(  I  m,  J2  T2  /  (N_ 
j=l 


*  i  o 

1 J 


- 1) 
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>  n*  T  T?  >  ao  n3 

—  L  i  —  *• 

ieSjj ■ 

where  0  <  a2  <  00  and  0  <  n*  <  00  is  a  constant  independent  of  the 

configuration 0  The  second  inequality  follows  from  the  fact  that  the 

quantity  under  the  first  square  brackets,  being  the  variance  of  a 

random  variable  which  takes  values  ci,C2,...,c  with  probabilities 

proportional  to  ,m,2 , 0 . . ,m. ^  is  bounded  below  by  a  positive 

constant  independent  of  the  configuration  and  i  e  S..f. 

J  J 

Thus  the  equations  (2.1.1)  and  (2.1.2)  hold  for  an  arbitrary 

k-1 

linear  combination  U.  =  E  c.R. . „  The  rest  of  Lemma  2.1  follows  from 

1  j=i  J  ^ 

Berry-Esseen  theorem  [3],  which  can  be  stated  as  follows. 


Theorem  B.  Let  Uj ,  U2 , 


0000 


be  independently  distributed  with  means 


E(Ui)9  E (U2 ) , o  o . .  and  let  b?,  3^  and  have  the  same  meaning  as  in 

n 

Lemma  2.1.  Let  F  (x)  denote  the  c.d.f.  of  E  {  U.  -  E(U.)}/S  ,  and 

n  1  i  n 

a_“  _L 

$(x)  the  standard  normal  c.d.f.  There  exists  a  constant  c  <  00  such 


that  for  all  x, 


n 


F  (x)  -  $(x) I  <  c  0  I  3.  /S3 
n  1  —  L  1  n 

i=l 


The  following  lemmas  have  been  proved  in  [lU]. 

Lemma  2.2.  (Mehra).  Let  V2  ,  .  .  .  .  ,  ),  n  =  1,2, - 

be  a  sequence  of  random  vectors  with  c.d.f.  F^n^(v)  (in  the  c-space  R^C^) 

depending  on  a  parameter  0,  and  let  =  f(V^n^)  =  E  c„V^n^  (f:  R^C^  ->  R^^) 

-  -  j=l  <3  J 

( n ) 

of  the  components,  such  that  as  n  -*  °°,  (i)  FA  ;(v)  converges  to  F  (v)  for 

each  0  (and  v),  where  FA(v)  is  the  c.d.f.  of  a  r.v.  V,  (ii)  for  every  linear 

-  -  -  0  -  - 

( n ) 

f,  GA  '(u)  converges  to  GA(u)  uniformly  in  0  (and  u)  where  GA(u)  is  the 
0  -  0  -  -  -  0  - 

c.d.f.  of  the  r.v.  f(V),  then  FAn\v)  also  converges  to  F  (v)  uniformly  in 

-  -  9  -  0  - 


I 


Lemma  2.3.  (Mehra)„  Let  the  sequences  {  }  and  {  F^^v)} 

n  =  1,2,,.,,  be  as  defined  in  Lemma  2 . 2  and  suppose  the  conclusion  of 
Lemma  2,2  holds  „  Let  f  i  R^C^  -*  R^  be  a  continuous  function,  then 


(n) 


.(n) 


(i)  the  distribution  G^  ;(u)  _of_  f(Vv  ')  converges  to  the  distribution 
G  (u)  of  f(V)  and  (ii)  the  convergence  is  uniform  in  0  ( and  u) 

Proof  of  the  Theorem  2.1.  It  follows  from  part  (i)  of  Lemma 
2„1  that  the  standardized  r.v.*s  Rls  R2 , »  •  »  .  ,R  in  the  limit,  as 
n  ->  °°,  have  a  (k-l)-  variate  normal  distribution.  This  verifies 
condition  (i)  of  Lemma  2.2.  Part  (ii)  of  Lemma  2.1  verifies  condition 

(ii)  of  Lemma  2.2.  Thus  the  convergence  of  the  standardized  Ri,R2,...,R 
to  the  normal  distribution  is  uniform  in  configuration.  Again,  the 


k-l 


statistic  W  ,  being  a  quadratic  form  in  Ri  ,R2  >•->»•  ,RV  is  a  continuous 
n  J£— _l 


mapping  from  R^  ^  to  R  ^ 


Thus  the  conditions  of  Lemma  2.3  hold. 


The  following  theorem,  proved  in  Rao  [17,  p«57]  shows  that  G.(u)  of 

O 

Lemma  2.3  is  the  X2  -  distribution  with  (k-l)  degrees  of  freedom.  The 
convergence,  by  Lemma  2.3  is  uniform  in  configuration. 

Theorem  C.  If  the  vector  R^=  (  Rj  -  E(R]_),  R2  -  E(R2 )  ,  .  .  .  „  ,R^  i  -  E(R^ 


has  a  normal  distribution  with  mean  vector  y  and  non-singular  variance- 


covariance  matrix  A,  then  the  quadratic  form  R’  A  1  R  has  a  X2,  _,(A) 

•vwv  ^  *vW  K— _L 

distribution  i.e.  a  non-central  chi-square  distribution  with  non-centrality 
parameter  A  and  (k-l)  degrees  of  freedom,  where  A  =  if  A  1  jn 

Remark  .  In  proving  theorem  2.1  we  assumed  the  aligned  observations 
to  be  distinct;  which  is  true  with  probability  one  if  the  aligned  obser¬ 
vations  have  a  continuous  c.dof.  In  practical  applications  of  the  test, 
equal  aligned  observations  can  be  allot ed  distinct  ranks  after  arranging 


them  in  some  order  at  random. 


r 
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An  alternative  method  of  ranking  will  be  to  allot  each  member  of 
sets  of  equal  aligned  observations  the  middle  rank.  For  this  method  of 
ranking  and  alignment  on  the  median,  the  conclusion  of  theorem  2.1  holds 
in  spite  of  the  fact  that  many  aligned  observations  may  become  zero.  The 
proof  of  Lemma  2.1  goes  through  so  long  as  there  are  distinct  positive 
and  negative  aligned  observations  in  each  block. 


CHAPTER  III 


ASYMPTOTIC  DISTRIBUTION  UNDER  THE  ALTERNATIVES 

3  o 1 .  Introduction 

In  the  previous  chapter,  the  asymptotic  null  distribution 

of  was  found  under  very  general  conditions  on  the  experimental 

layout  and  the  method  of  alignment.  In  this  chapter,  we  will  study 

the  unconditional  asymptotic  distribution  of  the  statistic  for 

certain  sequence  of  translation  alternatives  approaching  Hq  .  For 

complete  experimental  designs,  with  the  restriction  m. .  =  m.  for 

i  J  J 

all  i  (  for  each  j  ),  the  limit  distribution  was  studied  in 

Mehra  -  Sarangi  [lU].  The  results  of  this  section  extend  the  results 

of  [lU]  to  (i)  balanced  incomplete  block  designs,  when  m. .  =  0  or  1 , 

and  (ii)  partially  balanced  incomplete  block  designs  of  a  certain 

type  to  be  specified.  We  confine  our  attention  to  alignment  on  the 

mean  only.  The  unconditional  distribution  of  W  will  then  enable 

n 

us  to  compute  the  (Pitman)  efficiency  of  this  test  relative  to  the 
F  -  test,  the  competing  separate  ranking  tests,  and  certain  tests 
proposed  by  Lehmann  in  [ll]  and  [ 12 ] „ 

For  every  fixed  n,  consider  the  alternative  hypothesis: 

(3.1.1)  K  :  F.(  x  )  =  F(  x  +  0.  /  /n  )  (j  =  1,  2,...,k) 

^  J  J 

where  not  all  0^.?s  are  equal.  The  method  of  proof  is  the  same  as 
that  of  Andrews  [l]  and  depends  upon  certain  properties  of  U  -  statistics 
obtained  by  Hoeffding  [9].  But  for  the  sake  of  clarity,  the  two  cases 
mentioned  above  will  be  discussed  in  separate  sections. 


-  17  - 


3o2.  Balanced  Incomplete  Block  Designs 

Let  the  experimental  layout  be  a  balanced  incomplete  block 
design  with  parameters: 

n  =  number  of  blocks  in  the  layout , 
k  =  number  of  treatments  to  be  compared, 
m  =  block  size  i.e.  number  of  plots  in  each  block, 
r  =  number  of  times  each  treatment  has  been  replicated, 
p  =  number  of  blocks  in  which  any  two  specified  treatments 
occur  together. 

The  following  relations  are  well-known. 

(3.2.1)  n.m  =  r.k,  r(m-l)=p(k-l),  n  >_  k  >_  m, 

r  <  n 


As  before,  let  S.  denote  the  set  of  blocks  in  which  the 

J 

th 

j..  treatment  is  present,  and  S..,  the  set  of  blocks  in  which  both  the 

J  J 

and  the  treatments  are  present.  The  sets  S.  and  S. .  are 

J  J  J 

identified.  We  may  assume  (with  probability  one)  that 
(E)  All  aligned  observations  are  different  from  zero. 

The  conclusion  of  the  theorems  below  would  hold  even  without  this 
assumption.  We  are  making  such  an  assumption  specifically  to  overcome 
notational  difficulties. 


Let  Z. .  stand  for  the  aligned  observation  (under  assumption 
J 

E)  in  the  (i,j)^  cell  if  i  e  S  and  be  zero  if  i  i  S .  .  Let  r.  .  denote 

J  J  1 J 

the  rank  of  Z . .  in  a  combined  ranking  of  the  totality  of  aligned 

J 

observations  if  i  e  S .  and  r. .  be  zero  otherwise. 

J  ij 

If  Xi ,  X2,...o,X  be  m  independent  random  variables  distri¬ 
buted  identically  with  c.d.f.  F(x),  let  G(u)  denote  the  c.d.f.  of 


4 
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Zi  =  (  Xi  -  X  )  and  let  H(u,v)  denote  the  joint  Codcf.  of  (ZlsZ2)  where 


m 


X  =  (  E  X,  /  m).  Let  0.  denote  the  mean  of  the  0's  that  are  present 
i=l  1  1 

in  the  i..  block.  It  can  be  verified  by  direct  computation  or  from 
(l.Uol)  that  for  the  above-described  balanced  incomplete  block  design, 


(3.2.2) 


E(Rj )  =  l  r.  =  l 


E  (r. ,,/m) 


1J 


ieS 


j 


j?=l  ieS 


JJ' 


l  a..  =  Var  (R.)  =  T  (T?  /m  -  r?) 
N  JJ  J  L  i  i 

ieS . 

J 


°jj’  =  (W 


-  I  (T?/m  -  r?)/(m  -  l) 

ieS  1  1 

JJ 


_  k  9  k 

where  r.  =  I  r../m  and  Tf  =  E  r: 


=i  ij 


j=l 


i=l  ^ 


We  now  prove: 

Theorem  3.1.  If  the  hypothesis  be  true  and  the  following  conditions 


are  satisfied  ( i )  F ( x )  possesses  a  continuous  derivative  F  ’  ( x )  a. e . 

and  (ii)  there  exists  a  function  b(x)  such  that  [  {  F(x  +  0)  -  F(x)}  /0  ]  < 


b(x),  where  _/“  b(x)dF(x)  <  °°9  and  (iii)  for  every  m,  there  exists  a  function 

b  (u)  such  that  [  {  G^m^(u+  0)  -  G^m^(u)}  /  0  ]  =  [  {  G(u+  0)  -  G(u)}/0]  < 

m  -  — 

b  (u),  where  /  b  (u)dG(u)  <  then  the  statistic  defined  by  (1.4.2) 


m 


m 


converges  in  distribution,  as  n  -*  °°,  to  X2^_^  (A  )  where 


(3.2.3) 


A  = 
w 


3(m  -  1)/(1  -  3A)(k  -  1) 


{  G* (u) }2  du 


12  k 

I  ( e.-e )2 

J 

j=l 


—  k 

with  0  =  I  0./k  and, 
j=l  J 


A  = 


-oo  J  _ 


G(u)  G(v)  dH ( u , v ) 
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Let 


The  proof  of  the  theorem  follows  from  the  lemmas  proved  below. 


(3.2.4) 

(3.2.5) 


k 

=  /  3(k  -  1)  { -R  -  l  l  r  ,/m  }  /mn3/2  /(m_i) (i-3A ) 

J  j'=l  i  J 

=  -  /3(m  -  l)  {  [°°  (G’(u)}2  du  }  ( 0  -"0)/  /(k-l)(l-3A) 

'  —.00  J 


Further,  let  X  represent  the  variance-covariance  matrix 

||  6 .  ,  -  1/k  ||  with  =  1 ,2 , .  » .  ,k » 

i  J 

Lemma  3.1.  Under  the  assumptions  of  theorem  3.1,  the  random 

vector  T  =  (T  ,  .  T  .  )  converges  in  distribution,  as  n 

-  ^_n  nl  n2  3  nk  - 2 - - - 

to  the  k-variate  normal  vector  with  mean  vector  _y  =  (i+i,  yv)  anl 

— —  ■  -  -v-w  - 

variance-covariance  matrix  X 

Proof.  We  may  assume  without  loss  of  generality  that  the 
aligned  observations  are  all  different  from  zero  (see  assumption  E)„  It 
is  sufficient  to  prove  that  the  k2-statistics  given  by 


Tn( j j ' ) 


r . . /n3/2 
^  J 


(  =  1,2 


have  an  asymptotic  normal  distribution;  for  the  asymptotic  normality  of 
follows  by  a  linear  transformation,  The  later  result  is  derived  here  by 
arguments  parallel  to  those  of  Andrews  [l]  and  Mehra  -  Sarangi  [lH]. 


Define  6(y  ,y  ,)  =  1  if  y .  >  Yyi  Jy  Yy  ?  0  and 

.  k 

<S(yj,yj,)  =  0  otherwise.  Further,  define  hJ(y1,  y2,...,yk  )  =.?=1  ^(YyYy 


Let 


(3.2.6)  U 


,  ( j  j  '  ) 


n 

I  ^  ^i  ]_ 5  2 3 

V1  1  2 


n  n 

-  V  V 
L  L 


coo  /coo  oooo 


>Zi  k);  n 

V 


■ 
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It  is  easily  verified  that. 


The  R.HcS.  = 


n  n  n  k 

Ij  £  •  •  .  £ .  o  ,  .  .  *  o  £  £  5(Z  ,  Z  it)  /n 

V1  V1  V8JJ.  v1^1  J  J' 


:  V 

Ysjj* 


1  l  l  «(z,  .»  z  ,>)}  /  n2 


i.,=i  j>=i  j 


(3.2.7)  =  (  I 


rij  -  p}  /n2 


Define , 


f  i  1  ’  ) 

v J J  '  (z.  ,  z.  . . ,  z.  ) 

xl  x2  Xk 


( oi^  j  • «  t  aR ) 


h  (Za1l>  *  *  *  *  *  Z\k)  /k 


if  {  i1#  i2,...,iv  >n  S,.,  /  $ 


JJ 


=  0  otherwise 


where  (a  ,  a  , .  „  .  ,a )  extends  over  all  possible  permutations  of  the  set 

±  c.  K. 


of  integers  {  i  ,  i  i  }  such  that  a  e  S  $  being  the  empty  set 

1  2  k  J  JJ 

Let . 


U 


(jj')  _ 


^  ’  12  ’  0  °  5  ’'"k  ^ 


(Z.  ,  Z,  z  )  /[l 

k 


11  12 


where  (i  ,  i  ,o,.,i  )  extends  over  all  sets  of  integers  such  that 

jL  cL  K. 


1  -  il5  i2 * 


( i  i  * ) 

„ , o ,  i  <  n,  and  no  two  of  the  i*s  are  equal,  U  is  the 

K. 


average  of  (  ^  ),k!  terms  of  the  type  hJ  (Z  ,  Z  ,,,.,Z  ),  the  zero 

x  ot  X  cx_.2  a,  x 

12  k 

( ]  i  »  ) 

terms  of  which  however,  do  not  appear  in  U’  u  ,  It  is  easily  checked 


that  the  number  of  terms  of  the  type  h^(Z  ,,  Z  Z  ,  )  which  are 

a_  1  a_2’  a.  k 

12  k 
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( 1  i  * ) 

averaged  to  give  U  00  and  also  occur  in  U 


are  p{(k-l)!}.(  ) 


Let  ^  consisting  of  p[  n^  ^  j  ){(k-l)J}]  terms  denote  the  sum 


of  all  h^(Z  ,  Z  Q,..»,Z  )  terms  that  appear  in  the  expanded  form  of 

ot^-L  ot^k 


U 


(  i  i  »  )  (  J|  4  f  ) 

'  JlJ  hut  not  in  that  of  U  0  Then, 


U' 


(jj*)  =  I  nl  ck»  +  j(jj')  1/  nk 


k> 


It  can  now  he  verified,  as  in  [lU],  that 

Lim  {  Var  /n  (  -  u^j’)  )J  =  q 

n-x» 

/  •  •  y  \  /  •  •  y  \ 

Hence,  hy  the  following  lemma  of  Hoeffding,  if  either  /n  (  U  ^  ;-EU  ); 

* 

J,j'  =  1,2, ... ,k  or  /n  (U'^^  -  EU'^'^  );  j,j»  =  l,2,...,k  has  a 
limiting  joint  distribution,  the  other  has  the  same  limiting  joint 
distribution. 

Lemma  D.  Let  V^,  V  , ....  he  an  infinite  sequence  of  random  vectors 
V  =  (  V^\..„,  V^k^)  and  suppose  that  the  c.d.f.  F  (v)  of  V  tends 


n 


n 


n 


n 


n 


n 


to  a  c.d.f.  F(v),  as  n-*30.  Let  ^  ^  +  d^  ^  ,  where,  Lim  Ejd^^}2  =  0 

n  n  n  n-*a>  1  n  J 

for  all  j  =  l,2,...,k.  Then  the  distribution  function  of  Vf  =  (Vv  , .  . .  j'v  ) 

n  n  n 

/  •  «>  y  \  /  •  «  (  \ 

tends  to  F(v).  The  limiting  joint  distribution  of  /n  (  -  EU  ); 

j , j '  =  l,2,»..,k  is  proved  below  to  be  k2  -  variate  normal,  following  the 
theorem  of  Hoeffding  [9,  p.310]  on  U-statistics  for  non-identically  distri¬ 
buted  vectors,  which  can  be  extended  as  follows. 


Theorem  E.  Let  Z  ,  Z^, . . . ,2^  be  n  independent  random  vectors  of 
k  components,  having  the  c.d.f.  F_^(Z)  =  F  (  Z^\  Z^,,..,Z^k^ 


) .  Let 


22 


<j>^  (Z^,  Z2,cc„,Zh  ),  j  =  l,2,»0.,g  be  functions  syrnmetric  in  its  h  vector 
arguments  which  do  not  involve  n,  and  let 


J 


N  .  .  .  .  .  (Z.  )  =  E  {  <J>J(Z.  ,.,,,Z.  )  ( 

l(i  )  .1  1  1  l  ’  ’  l,  '  if 

v  12  v-1  v+1  h  v  1  h  v  J 


I 


I 


1 


-  E  {  <*>J(Z  ,...,Z  ) 

1  h 


^(v)(  zv  >  = 


(Z  )  (v  =  1,2, . . ,n) 

h-1 


where  ( i^ ,  i^ , „ . „ , i^  ^ )  extends  over  1  <_  i  ,  i^,...,  i^  ^  n ,  with  no 
two  i's  equal  and  no  i  equal  to  v .  Suppose  there  exists  a  number  A 
such  that  for  every  n=l,2,...  and  j=l ,2 , . . . ,g, 


[  {  <|>J  (Z1#Z2,...,Zh)  }2  dF.  (Z1)90sdF.  (Zh)  <  A 

1  h 


(  1  <_  i  ,  i2 , .  . .  ,  i^  <_  n ,  no  two  i  '  s  equal ) 


E  |  ’1'J1(V)(ZV)  |3<  “  (v  =  1,2, ... ,n;  j  =  1,2,. ..,g)  and 

3/2 


n  g 


Lim  l  E|  I  c.-r5.,  x(Z  )  |3/ 
n-*»  v=l  Joi  J  ^  v  / 


XE  1JLCJ?Jl(v)(Zv) 


=  o 


for  arbitrary  constants  (c_,  c_,.„.,  c  ).  Then,  as  n-*»,  the  joint  dis- 

12  g 

•  •  « 

tribution  of  { (  -  EU^  )/  /Var  },  j  =  l,2,..„,g  tends  to  a  g- variate 

normal  distribution  with  null  mean  vector. 

If  h  =  k,  and  g  =  k2 ,  it  is  easily  verified  in  our  special  case 

that , 


J...J  {  <tfJJ,)(z1,z2,...,zk)  p  dF.  (z1)  dr.  (z2)...dF,  (zk) 


(k-l) 


. 
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E 


( jj  ’  ) 
l(v 


3  <  8 (k-l) 3 


oo 


It  follows  from  above  that, 


n 


(3o208)  l  E 


v=l  (j,Jf) 


y  c  ) 

jj  ’  Mv)  (v 


3  <  M  n 


where  0  <  M_  <  °°  is  independent  of  n.  Consider  v  e  S .  .  , „  If 
1  JJ ' 


Z-  »Z .  900c ,  Z  .  are  (k-l)  arbitrary  vectors  different  from  Z  ,  then 
1  2  1k-l  v 


(jj  '  ) 


riZu  i  i  <V  =  E  {  (Z  z  ...z  )  I  z  } 

MVA112">  lt-1  1  Xk-1  V 


-  E  {  *(JJ,)  (Z  Z ...... Z  )  f 

1  k-l 


k  r 

=  l 

j'=i 


I  1 E(  s( 

i.S7^v  J  J 

J 


l 


VueSj,r 


f  E(6(  Zi.j’Zvj’)lZv)'E6(Zi.j’Zvj')  f 


/k! 


Thus, 


?i(J5>  (V 


k 

Y 

L 

'=1 


v 

L 


(  E({(  Zvj’ZiJ|)|Zv) 


E6  (ZvJ,Z.JT)  } 


I  {E(6(Z.J,Zv.,)|Zv)-E6(Z..>Zv.,)  } 

jj  ' 


/k! 


rk  {  G(Z . . )  -  1/2  }  +(r  or  p)  [  {  1/2  -  G(Z  ,)}  +  0(/n) 

j*=!  VJ 


=  n 


/  k: 


-1/2, 


m  {  G(Z..)  -  1/2  }  +  (m/k  or  m(m-l) /k(k-l )  )  T  {  l/2-G(Z  ,)  +  0(n  )/k* 

ij  j '=!  VJ 


where  the  first  or  the  second  alternative  is  chosen  according  as  j  =  j  * ,  or 


-  2k 


j  ^  j ’ .  It  is  seen  from  above  that  the  variance  of  any  linear  combination 
of  the  left-hand  side  is  bounded  away  from  zero.  Thus,  we  can  find  M0  >  0, 
such  that , 


n 


J’U,,  *3AN  (v  }2-v 

v  =1  (j  ,j  '  ) 


which  with  (3.2.8)  verifies  the  last  condition  of  theorem  E.  Thus  all 

( 1 1 ?  )  ( 1 1 '  ) 

conditions  of  the  theorem  are  satisfied,  proving  that  {  U  -  EU  dd  }/ 

_ _  /  •  •  f  \ 

/  Var  U  ^  ;  j , j  T  =  l,2,,..,k  have  a  joint  k2  -  variate  normal  distribution, 

as  n-*».  The  asymptotic  joint  normality  of  {  UJ'^  ^  -  EU *  ^ ^  ^}//  Var  U'^  ^ ; 
j , j  *  =  l,2,...,k  follows  from  the  results  on  p.21,  using  the  fact  that 

,  ( J  j  ’  ) 


n.Var  U1 


exists  -  which  is  verified  below.  The  mean  vector  and  variance- 


covariance  matrix  of  the  limiting  normal  distribution,  following  the  results  of 
Hoeffding,  can  now  be  calculated  by  computing  their  values  for  a  fixed 
n;  and  then  proceeding  to  the  limits  as  n  tends  to  infinity.  The  existence 
of  the  mean  vector  and  the  variance-covariance  matrix  and  their  limiting 
values  have  been  shown  below. 

n  k  . 

Using  the  relation  r. .  =  I  I  6(Z. . ,  Z  )  +1,  it  is  easily 

1J  U=1  V=1  Ij  UV 

verified  that, 


E(r . . )  =  nm/2  -  /n.m 
1  J 


{  G’(u)}2  du 


(  ej  -  ei  )  +  0(1) 


Thus, 


l  E(r.  .)  =  nmr/2  -  >/npk 
1  *1 


ieS .  .  , 
JJ 


°°  {g’(u)}2  du 


J  —  oo 


(  0J  -  0  )  +  o(n) 


£  £  E(r .  .  ,  )=nmp  (k-l)/2  +  /npk 


ieS  ,  jVj=l 


{ G 1 (u) }2du 


(0  -0)+O(n) 


Using  (3.2.L)  and  (3.2.5),  this  gives 


(3.2.9)  Lim  E  (T  .)  =  y.  =  -  /3(m-l) 
ft-KC  ^  J  <J 


{G ’ (u)  }2du 


(  0  -0)/ Ak-l)  (1-3A) 


Similarly, 


n.  k 


n  k 


Var  (r.j)  =  E  {  £  £  «(  Z.,,  Zj  } 2  -  E2  {  l  l  6(Z.,,Z_)  } 


L  ~  '  i  j  5  uv‘ 
u—l  v=l 


U=1  V=1 


1J  uv 


V  V  V  y 
L  L  L  L 


1  '  U  S  V 


E{  6(Z  ,Z. ,  )6(Z.  .  ,Z  )  ) 

1  i j  *  i’s  i.i  ’  uv  J 


ij 


E{  6(Z.  .,Z.(  )  }E{  6 (Z  ,Z  ) J 

1  IJ  1  '  S  J  1  lj  uv  1 


=  n2m2/12  +  0(n^^) 


Cov  (r  ,  r..,)  =  n2m2(A  -  1/U)  +  0(n~^2) 
ij 


Thus , 


Lim  (  Var  T  . ) 

_  nJ 

n-*» 


=  Lim  {  3(k-l)  }/{  n3m4 (m-1 ) (l-3A ) } .Var  { (m-1 )  \  r. 


n-x» 


ieS 


ij 


k 

V  V 

L  L 


jV j=l  ieS 


r  , } 
ij 


JJ 


=  Lim  {3(k-l)/n3mf(m-l)  (l-3A) } .  {m(m-l  )r  (m2n2/12  )-m(m-l  )r{m2n2  (  A-l/4 
n-*» 


(3.2.10) 


=  (l-l/k) 


Since  the  expansion  of  Var{  (m-l)  \  r_  \  \  r  _  ,  } 

ieS .  j  Vj=l  ieS  , 

J  J  J 


using  the  formula  Var  (EX.  )  =  E  Var(X. )  +  2  EE  Cov(X.,  X. t )  gives 


i=l  x  i=l 


i>i 


f  11 


m(m-l)r  terms  of  each  of  the  type  Var(r  )  and  Cov(r.s,  r..,)°  The  only 


ij  ij  ij 

other  terms  of  the  type  Cov(r^.,  r_^ ,  .  f  )  when  i  ^  i',  cancel  out.  The 
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limiting  covariance  term  is  similarly  calculated  or  follows  from  the 


k 


equation  E  T  ,  =  0  coupled  with  the  fact  that  the  above  limit  is 
j=l  nj 

independent  of  (  0^,  02,„..,  @k  ),  so  that, 

(3.2,11)  Lim  Cov  (T  ,  T  f)  •«  -  Lim  Var  (T  )  /  (k-l)  =  -1/k 

n-*» 


A  proof  of  asymptotic  normality  together  with  results  (3.2,9),  (3.2,10), 

and  (3.2,11)  completes  the  proof  of  Lemma  3.1. 

Lemma  3.2.  Under  the  assumptions  of  Theorem  3.1,  the  conditional 

variance-covariance  matrix  of  T  =  (  T  n ,  T  _ , , , , ,T  ,  ),  given  by  (3.2,2) 

- — - —  n  nl  n2  nk  - - - 

'W 

tends  in  probability,  as  n-*»,  to  the  unconditional  variance-covariance 
matrix  given  by  Lemma  3.1. 

Proof.  Since  the  expectation  of  the  conditional  variance-covariance 
matrix  is  the  unconditional  variance-covariance  matrix;  the  conclusion 
of  the  lemma  will  follow  from  Chebyshev’s  inequality  if  we  prove  that  the 
variance  of  each  element  of  the  conditional  variance-covariance  matrix 
tends  to  zero,  under  the  alternatives,  as  n-*»,  This  is  proved  as  follows. 


{  l  r_  /n3/2  }  =  l  {l  r^/m  -(  \  r_/m)2  }/n3 


ieS 


Var  d .  .  =  Var 
JJ 


<  2 


ieS  j=l 
J 

k  k 

l  {  l  r?j/m  -  (  l  r  /m)2  } 

ieS  j=l  j=l 

J 


j=l 

/n6 


Var  (  J  l  r?  }/m2  +  Var  {  l  (  l  r  )2}/m4 

ieS,  j=l  J  ieS  j=l  J 

J 


J 


/n 


-*  0,  as  n-x=°, 


since  each  of  the  two  terms  above  tend  to  zero  by  virtue  of  the  properties 
(i)  Cov  (  r2k,  r2k  )  =  0(n3)  ,  (ii)  Cov(  r2^  ^2ir22  ^  =  » 


' 


t 


— . -  ■ 
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( iii )  Coy  (  r^r^,  ?  2^r  22  )  -  0(n3) 


Proof  of  Theorem  3.1.  Let  T*  =  (T  ,  .  T  „ . .  . .  «T  ,,  ,  s  )  and 
- -  n  nl’  n2*  ’  n(k-lj 

/w 

A*  =  ||  <5  .  .  f-  1/k  ||;  j , j f  =  1,2,,,,, (k-l)  denote  their  associated 
^  J  J 

variance-covariance  matrix.  Then  by  virtue  of  Lemma  3.2,  using  the 
following  theorem  given  in  Cramer;  both  W  and  T*  T* '  have  the 


n 


n  ** 


n 


same  asymptotic  distribution.  Following  theorem  C  [ IT ]  the  latter  has 
a  non-central  chi-square  distribution  with  (k-l)  -  degrees  of  freedom 

i  k  9 

and  non-centrality  parameter  A  =  y*  A*”-1  y*  ’  =  I  y.,  where 

V'V*/VV'  AM  9  1  1 

j=l  J 

y*  =  (  y  ,  y  , .  .  .  ,  y  )  -  which  agrees  with  the  value  given  by  (3.2.3). 
1  d  K-l 

This  proves  the  theorem. 

Theorem  F.  Let  C  ,  be  a  sequence  of  random  variables,  with 


the  distribution  functions  F  ,  F^, 


e  ©  9  o 


Suppose  that  F  (x)  tends  to 


F(x),  as  n-*».  Let  r)  ,  ^2 , . . . .  be  another  sequence  of  random  variables, 


and  suppose  that  nn  converges  in  probability  to  a  constant  c.  Putting 

Z  =(c  /  n  ),  if  c  >  0,  the  distribution  function  of  Z  tends  to 

n  n  n  n 

F(cx) . 

3.3.  Partially  Balanced  Incomplete  Block  Designs, 

The  results  of  this  section,  contained  in  Theorem  3.1(A)  below, 

are  direct  extensions  of  the  results  of  Theorem  3.1.  Let  our  experimental 

layout  be  a  partially  balanced  incomplete  block  design  with  parameters 

(n,  k,  m,  r,  p..,)  where 
J  J 

p  .  .  ,  =  number  of  blocks  in  which  both  the  and  j 

JJ  ’ 

treatments  occur  together, 

and  the  other  symbols  have  the  same  meaning  as  in  section  3.2. 
We  further  assume  that. 


Lim  t  Pjj  »  /nf  =  Pjj  ?  >  ° 


(3.3.1) 
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Then  the  equations  (3.2.1)  can  be  modified  as. 


(3.3.2)  n.m  =  r.k,  r(  m  -  1  )  =  £  p  ,  ,  n>_k>_m 

J  J 


J 


r  <  n 


We  now  state, 

Theorem  3.1(A).  Assume  for  each  n,  the  truth  of  the  hypothesis 

K  ,  and  that  the  conditions  of  Theorem  3 . 1  and  the  condition  (3.3.1) 

are  satisfied.  Then,  for  the  above  partially  balanced  incomplete  block 

design,  the  statistic  W^  defined  by  (1.4.2),  converges  in  distribution, 

as  n-*30,  to  X2.  .  (  A*  ),  where, 

—  —  k-1  w  - 


(3.3.3) 


3(m-l) (k-1 ) /k2 ( 1-3A 


) } 

r 

’  G'(u) 

2  i 

}  du 

—00 

— 

2 


x 


0 

J 


k 

k/m(m-l)  £ 

j  ’/j=l 


2 

} 


Proof  Let 


(3.3.4)  Tnj  =(v/3(k-l)/m/(m-l)(l-3X)  } 


k 


R,-i/m  I  l 


1  ^  ^  rij’ 

J  j’=l  ies.  .  ,  J 

JJ 


>  /n 


3/2 


and 


(3.3.5)  y .  =  -f/3(m-  l)(k-l)/k/l-3A  } 

J 


2  -I 


r 


{  G’(u)}  du  / 0 .  -  k/m(m-l) 

J 


.  I 

3  '/j=l 


where  S . . , 
jj  ' 

Lemma  3.1, 


is  the  same  as  in  Lemma  2.1.  Following  the  lines  of  proof  of 
we  can  prove  that  the  random  vector  T^  =  (T  ,  T .  .  .  ,T^^) 


has  an  asymptotic  k  -  variate  normal  distribution  with  mean  vector 


_ 

V 


-  29  - 


y  = (ui SU2 > ° * • »VU  )  and  variance  -  covariance  matrix  A.  It  can  be  proved  as 

in  Lemma  3.2  that  the  conditional  variance  -  covariance  matrix  of  T  ,  given 

_ n 

by  (3.2.2)  tends  in  probability,  as  n-*»,  to  the  unconditional  variance- 
covariance  matrix  JV,  The  proof  of  the  theorem  thus  can  be  accomplished 
exactly  in  the  same  manner  as  Theorem  3.1  and  need  not  be  repeated  here. 

3.4  "Balanced"  Incomplete  Blocks  With  p  Observations  Per  Cell 

The  results  of  this  section  are  an  extension  of  the  corres¬ 
ponding  results  of  section  3.2.  Suppose,  in  the  experimental  design  of 
section  3.2,  p  observations  are  available  for  every  treatment  that  is 

present  in  a  block  i.e.  m. .  =  o  or  p.  The  notations  of  this  section  have 

ij 

the  same  meaning  as  the  corresponding  notations  of  section  3.2.  Further, 

let  D. .  denote  the  sum  of  the  ranks  under  the  treatment  in  the  i"^ 

ij 

block.  Define 


(3.U.1) 


nj 


/ 3  ( kp-l)/fm ( np )  ^ ^ 2  /(mp-1 )  ( 1-3A  )  } 


R 


k 

-  I 

j 


'=1 


and 


(3.4.2) 


y  =  -  (  (m-l)  /3p(kp-l) 
J 


} 


{  (k-l)  /(mp-l)(l-3A)  } 


X 


r°°  2 

{  G '  ( u )  }  du 

. 

_  — oo 


for  j  =  l,2,...,k.  It  is  easily  verified  that 
equivalent  to  the  statistic  of  (1.4.2). 

Theorem  3.1  (B).  Assume  for  each  n  the  truth  of  the  hypothesis  K  , 

and  that  the  conditions  of  Theorem  3d  are  satisfied.  Then  for  the 

above  experimental  design,  the  statistic  defined  by  (1.4.2) ,  converges 

2 

in  distribution,  as  n^° ,  to_  x ^  p )  »  where 


V 


V  2 

l  T  • 

>1  nJ 


is  asymptotically 


(3.U.3) 


£  2 
4w=I  »J 

J=1  J 


2  /  2 
=  {  3p(kp-l) (m-l)  }  j  {  (mp-l) ( 1— 3X ) (k-l)  } 


{  G'(u)  }  du 


I  <6j 
J=1  J 


-  0) 


Proof,  The  asymptotic  normality  of  T  =  (T  „ ,T  ,  ),  its 

-  n  nl  n2  nk 

limiting  mean  vector  and  variance  -  covariance  matrix  are  deduced  from 
the  results  of  section  3,3c 

This  is  done  by  associating  separate  0.’s  with  the  p  observations  in  a 

J 

cell  and  noting  that  our  experimental  design  reduces  to  a  partially 

balanced  incomplete  block  with  two  associate  classes,  having  parameters 

(n,  kp,  mp,  r,  p  )  where 

J  J 


n 


=  number  of  blocks 


kp  =  number  of  treatments 


mp  =  block  size 


r  =  number  of  replications  of  each  treatment 

p..,  =  number  of  blocks  in  which  both  the  j^h  and  j 
J  J 

treatments  occur  together. 

This  is  r  or  p  according  as  j  and  j*  do  or  do  not  occur  in  the  same  cell, 

The  joint  distribution  of  T  =  (T  _ ,T  _ , . „ „ ,T  .  )  follows  by  a  linear 

n  nl  n2  nk 

transformation  and  the  results  of  Theorem  3.1(B)  follows  as  before  after 
setting  the  0  * s  within  the  same  cell  to  be  equal. 

Remark  It  may  be  observed  that  in  Theorems  3.1,  3.1(A)  and  3.1(B), 


the  asymptotic  distribution  of  the  statistic  W  is  the  same  as  the 

ft  2\ 

asymptotic  distribution  of  £  T  „  The  later  expression  with  A  replaced 

\j=l  nJ  1 

by  a  consistent  estimator  may  therefore  be  proposed  as  an  alternative  test 


statistic.  Many  such  consistent  estimators  of  A  are  provided  by  the  result 
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of  Lemma  3.2  (or  its  analogues  in  Theorem  3.1(A)  or  3.1(B)),  This  statistic 
may  more  easily  be  applied.  However,  whereas  the  convergence  of  to  the 
chi-square  distribution  is  uniform  in  configuration;  the  same  property  may 
not  hold  in  case  of  the  later  statistic. 


CHAPTER  IV 


COMPARISON  OF  EFFICIENCIES 


^ .1  The  Concept  of  Asymptotic  Efficiency 

The  concept  of  asymptotic  (local)  efficiency,  originally 
given  by  Pitman  and  subsequently ■ studied  by  Noether  [15],  seems  a 
satisfactory  criterion  for  comparing  two  consistent  test  procedures . 
Suppose  S  and  S*  be  two  consistent  test  procedures  and  {k^}  be  a  sequenc 
of  alternatives  approaching  the  null  hypothesise  For  a  fixed  n  and  a 
pre-assigned  level  of  significance  a,  let  N^(n)  and  N2(n)  respectively 
denote  the  sample  sizes,  required  to  give  the  same  power  3  to  both  the 
tests  S  and  S* .  The  asymptotic  efficiency  of  S  relative  to  S* ,  for  the 
sequence  of  alternatives  {k^}  ,  is  defined  to  be 

®s,s*  =  es,s*  (“’6’  {kn})  =  Lim  N2(n>  7  Nl(n) 

’  9  n-**> 

if  the  limit  of  the  right-hand  side  exists.  In  practice,  the  limit  is 
often  independent  of  the  choice  of  (a, 3);  but  the  sequence  of  alterna¬ 
tives  have  to  be  suitably  chosen  to  find  its  value . 

^ 0 ^  Evaluation  of  Asymptotic  Efficiency 

In  this  section  we  give  exact  expressions  for  the  asymptotic 
(Pitman)  efficiencies  of  our  tests  relative  to  the  competing  test  pro¬ 
cedures,  namely,  the  classical  F-test  and  the  Durbin's  (or  appropriate 
separate  ranking)  test  procedure.  There  will  be  some  overlap  of  the 
results  stated  here  with  those  in  section  h  of  Mehra  and  Sarangi  [ lU ] 0 
It  has  been  used  by  Andrews  [l]  and  pointed  out  by  Hannan  [6]  that  if 
under  the  alternatives  both  the  test  statistics  have  non-central  chi- 
square  distributions  with  the  same  degrees  of  freedom,  then  the 
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asymptotic  relative  efficiency  of  the  corresponding  test  procedures  is 
given  by  the  ratio  of  the  corresponding  non-centrality  parameters.  We 
discuss  each  of  the  previous  cases  separately. 

I .  Balanced  Incomplete  Block  Design 

The  classical  F-statistic  is  distributed  in  the  limit  as  n-*», 
as  Xk-1^AF^  where  , 

k 

(4.2.1)  A  =  (m-1)  {  l  (6.  -  0)2  }  /  (k-l)o2 

j=l  J 

if  the  variance  a2  for  the  parent  distribution  exists.  The  corres¬ 
ponding  Durbin's  test  statistic  W^*,  based  on  independent  ranking  within 
each  block  is  given  by 

k 

W  *  =  { 1 2/kp (m+1 )  }  l  {R .*  -  £  (m+1)}2 

n  j=l  J 

Its  large  sample  behaviour,  as  n-*»,  is  given  by  the  following  theorem. 

Theorem  G.  Assuming,  for  each  fixed  n,  the  truth  of  the  hypothesis 
K^,  if  (i)  F(x)  possesses  a  continuous  derivative  F'(x)a.e.  and  (ii)  there 
exists  a  function  b(x)  such  that  [{F(x+0)  -  F(x)}  /  0]  <_ b(x)  where 

•CO 

b(x)dF(x)  <  then  the  statistic  W  *  converges,  as  n-*=°,  to  y}  ( A  „), 

n  k-1  w* 

*  ~oo 

where 

k 

(4.2.2)  A  =  {12m(m-l)/(m+1)(k_1)}  {  [F9(x)]2}2T  (0,  -0)2 

w  J-»  j=i  J 

The  results  of  (4.2,1)  and  the  preceding  theorem  are  contained  in 
van  Elteren  and  Noether  [19]. 

From  (3o2„3),  (4,2.1)  and  (4.2.2),  taking  the  ratio  of  the 
non-centrality  parameters,  we  get 
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(4,2.3) 


w,w 


* 


(m+l) 


/4m(l-3A) 


r00 

{ G f (x) }2dx  / 

rco 

{F* (x)}2dx 

• 

• 

-00 

00 

(4.2.4) 


/(1-3X) 


r°° 

{G* (x) }2  dx 

— oo 


2 


where  A  is  defined  by  (3.2.3). 

Explicit  evaluation  of  the  expressions  (4,2.3)  and  (4.2,4) 
seems  difficult  in  most  cases;  except  for  the  normal  form  of  the  parent 
distribution,  for  which  we  give  computations  below.  For  the  special 
case  where  m=2,  it  is  easily  checked  that  A  =  1/6  and  the  above 
expressions  can  be  evaluated.  Some  values  of  e^  p,  assuming  different 
forms  of  F  are  given  bdiow. 


F 


Normal 

.955 


Uniform 

.889 


Exponential 


Double 

Exponential 


1.5  1.172 


The  values  of  e  r  from  the  above  table  indicate  the  rather  high 

w,F 

efficiency  of  our  test  relative  to  the  F-test. 

Efficiency  for  Normal  Parent  Distribution 

Suppose  F’(x)  =  (/2ttc0  ^  exp  {  -  tj(X  -■■-■) 2  } 

Then,  by  direct  computation,  we  get 

G*(x)  =  {/m/a  /27T(m-I)}  exp{-  mx2/2(m-l)a2 } 

f  oo 

{F'(x)}2dx  =  and 

j  -00 


oo 


—  oo 


{G*(x)}2dx  =  /m"  /  2a/Tr(m  -l) 
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If  x  (i=l,2,3;  j=l ,2 , , , „ ,m)  be  3m  independent  random  variables 

J  _ 

_  n  m 

distributed  identically  with  c.d.f.  F(x)  and  x.  =  —  T  x. then 

1  m  j=l 


A  =  G(u)G(:v)dH(u,v)  =  Pr[  (xn-  x:  )  -  (x21  -  x2 )  >  0, 

J  « 

—  00  —  CO 

(x12-  Xi)  -  ( x 3 ! —  x 3 )  >  0] 

=Pr[u>0,v>0] 

where  (u,vj  have  a  joint  distribution  that  is  bivariate  normal  with 
zero  means,  common  variance  {2a2(m-l)/m}  and  covariance  {-a2/m}.  Thus, 


A  = 


m 


0  o  2ttct2  A(m-l)2-l 


m(rn-l)  / 

exp<  -  - 

c2A(m-l)2-l]  V 


x2+  y2+  — M — ' 
(m  -1), 


dxdy 


m 


tt/2 


2tt  A(m-1  j^-1  0 

tt/2 


exp 


m(m-l )r2 
/  _  _ 

4(m-l)2-l 


=  A  (m-l )  z-1/^tt 


=  /+( m-l)2-l/2 


TT 


d0 


0 

tt/2 


(m-l)+Sin0  Cos0 


d0 


2(m-l)+Sing 


.  Sin0  Cos  0, .  , 

1  +  (m-l)  ,  >rdr 


de 


—  arctan  /  ( 2m-3 ) / ( 2m- 1 ) 

TT 


Substituting  these  in  (4.2.3)  and  A  .2  <.4),  we  get 


(4.2.5)  J 


w,w*  (normal)  =  (m+l)/j+(m_i)J  1—  (  ^ /tt )  arctan  /( 2m-3 ) /( 2m- 1 ) 


-1 


w 


, F  ( normal )  = 


3m 


Arr(m-1 ) 


1—  (  3/tt  )  arctan  A  2m- 3 )  /( 2m-l) 


-1 
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It  is  easily  checked  that  (4.2.5)  can  be  rewritten  as 


[  eW)W*  (normal)  =  ^  3a  V^(l+a2) 


(4.2.6) 


e  r  ( normal )  = 

w,  r 


3(3-“2)A.(i+«2) 


l-(  /tt  )  arctana 


1-(3/tt)  arctana 


-1 


-1 


where  a2  =  (2m-3)/(2m-l) .  Differentiating  by  the  product  rule, 


d_ 

da 


{ e  „ ( normal ) } 
w,w* 


12aTT/{  (l+a2  )  (  tt— 3  arctana)}2  I  (5-3a2)/l6  + 


+  arctan  a- 


since  the  quantity  within  the  square  brackets  has  a  negative  derivative 

with  respect  to  a  and  thus  has  a  maximum  when  a  has  the  minimal  value  of 

^ / /3  (for  m=2)  -  which  can  be  verified  to  be  negative.  Thus  e  „( normal) 

w,w* 

given  by  (4.2.5),  attains  a  maximum  of  3/2  for  m=2,  and  monotonically 
decreases  to  1,  as  m  tends  to  infinity.  Similarly, 


d_ 

da 


( normal ) } 


8 


9/ {  4a  ( l+a2  ) 2  (  tt-3  arctana)2}  (3-a2)/a-  U7T/3  + 


+  8  arctana"] 


The  quantity  within  the  square  brackets  can  be  verified  to  be  positive  for 
m=2  when  a  =  "V/3,  negative  for  m=3  when  a  =  /3/5  and  its  derivative  with 
respect  to  a  can  be  checked  to  be  negative.  Thus,  considering  only  integral 
values  of  m,  e  .-( normal)  given  by  (4.2.5)  attains  a  minimum  of  3/tt  for  m=2, 

W  ,  r 

a  maximum  for  m=3  and  monotonically  decreases  to  3/ tt  as  m  tends  to  infinity. 
Thus, 


e  r  >  3 / tt  for  all  m. 

w,F  - 


(4.2.7) 
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The  following  table  gives  some  values  of  the  efficiency  expressions  (4.2,5) 


m 


3 

1°  355 

e  r( normal)  3/tt=  0 95^+9  .9662 

w ,  r 


e  u,  ( normal )  1.5 

w,w* 


4 

1.263 

,9646 


5 

1 0  210 

.9631 


1 

3/tt 


The  behaviour  of  e  #( normal)  verifies  the  intuitively  clear 
fact  that  the  gain  in  efficiency  due  to  inter-block  comparisons,  by  using 
the  joint  ranking  procedure  instead  of  the  separate  ranking  one,  decreases 
to  zero  as  the  size  of  the  blocks  tends  to  grow  larger,  A  statement  in 
these  lines  was  made  by  Hodges-Lehmann  ([8],  Pp  495). 

Efficiency  for  large  m 


The  preceding  remark  that  e  „  tends  to  1,  as  m-*»,  if  the 

w,w* 

parent  distribution  is  normal  holds,  in  fact  for  all  cases  where  the 
parent  distribution  has  a  finite  expectation.  This  result  was  proved  in 
[l4]  and  we  give  a  proof  here  for  convenience. 


Theorem  4.1  If  (i)  for  the  parent  distribution  F(x), 


xdF(x)=a  <  00 


and  (ii)  there  exists  a  mn  such  that  for  all  m  >  mn,  there  exists  a 


function  b(u)  such  that  [{Gm(u+0)  -  Gm(u)}/0]  =  [ (G(u+0 )-G(u) }/q ] <  b(u), 


m. 


where 


b2(u)dti  <  00  ,  then  e  #(F)  given  by  (4,2.3)  satisfies 


w,w 


00 


Lim  e  „(F)  =  1 
w,w* 


n-*°° 


Proof:  Suppose  x„.  (i  =  1,2,...;  j  =  l,2,...,m)  be  independent  r.v.rs 

J 


1 


m 


with  identical  c.d.f,  F(x)  and  let  x.  =  —  Ex...  Since  F(x)  has  a  finite 

1  m  j=l  1J 

expectation,  x_^  tends  to  a  almost  surely  for  all  i  by  the  strong  law  of 
large  numbers  and  so 

Lim  A  =  Lim  C(xn-  x})  >  (x21~  x2 )  ,  (x12-  x1)>(x31-  x3 )  ] 
m^°°  m^°° 


=  Pr[xn  >  x21,  x12  >  x31]  =  1/4 


' 
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Further, 


Lim  G(x)  =  Lim  Prtx^-  x1  <_  x]  =  F(x  +  a) 
m-*30  m-+°° 


which  implies,  using  condition  (ii)  and  the  dominated  convergence  theorem, 

f  00  #oo 


Lim  [  {G’(x)}2dx  /  {F’(x)}2dx]  =  1 


m-*» 


00 


00 


This  proves  the  theorem. 

Remark;  The  following  example  shows  that  the  condition  of  finite¬ 
ness  of  the  expectation  is  essential  for  the  conclusion  of  Theorem  4.1. 
Suppose,  the  r.v,  has  a  Cauchy  distribution  of  the  form 


F' (a,b,x)  =  — 


b2+(x-a) 2 


,  b  >  o 


so  that  the  expectation  does  not  exist.  Since  it  is  easily  seen  that 

e  *  as  given  by  (4,2,3)  is  independent  of  any  change  of  scale  or 
w ,  ww 

location  of  the  r.v,,  we  assume  without  any  loss  of  generality  that  a=o 
and  b=l.  In  this  case,  (x^-Xj)  is  a  Cauchy  variable  having  density 
function  F'  (0,2(m-l)/m,x)  which  tends  to  F'(0,2,x)  as  m-*30.  So 

Lim  G* (x)  =  F ' (0 ,2 ,x ) . 
nr*30 


It  is  easily  verified  that, 


Lim 

m-*» 


{  G  *  ( x ) }  2  dx  / 


{  F T  (x)}2dx 


=  1/2 


Lim  A  =  Lim  Pr  [(x^-  x^ )  >  (X21-  x2 )  ,  (xi2-  xi  )  >  (x31-  X3)J 
m-*00  m-*30 

m  1  m  1  Xl2 

=  Lim  Pr  [  (  £  x  . )/m  >  (l-— )x2i~  (  ^  x  )/m  -  (1-— )xi 


m-^°° 


j  =  3  1J 


L  "21 

j=2  J 


m 


m 


(  ,Ix1j)/m>  (1-J)x3l-  <1  x3j)/m-  (1-^)xi2+ 

J  3  I  2 


m 


1 


X11 


=  Pr  [  x  >  u,  x  >  v] 


' 


t 
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where  x,  u,  v  are  independent  Cauchy  r.v.'s  of  the  type  F* (0,1 ,x) ,F' (0 ,3,x) 
and  F’(0,3,x)  respectively.  The  right-hand  side  can  he  simplified  to 


r  oo 


1/4  +  1/tt^ 


{(arctan  ^-)  2/(  1+x2  )  }dx 


Substituting  in  (402=3),  we  get 


r  oo 


(4„2„8)  Lim  e  „  (Cauchy)  =  [4 

w,w* 

m-x» 


48 

3  J_ 


IT 


{(arctan  ^-) 2/ ( 1+x2  )}  dx]  1 


For  evaluating  the  integral,  consider 
1(3)  = 


{(arctan  —f/  (l+x2)}  dx  ,  3  >  1 
•a  oo  p 


Let 


KB)  =  2 


{x  arctan  (-x/3)}/  {(l+x2)(32+  x2)}dx 


•'—00 


To  find  l(3),  consider 


L(B)  “ 


{z  log  ( 3-iz ) }/{ (l+z2 ) ( 32+z2 ) )dz 


T 


where  refers  to  the  contour  bounded  by  the  x-axis  and  the  upper  half 
of  the  circle  with  centre  at  the  origin  and  radius  R.  Calculating  the 
residues  at  z  =  i,3i;  since  the  integral  vanishes  along  the  circular  arc 
as  we  get 


f  00 


x  log(3-ix)  ,  Itt  3+1 

-  dx  =  -  log 

2  W  o2  .  „2  N  r>  9  n  2P 


( 1+x2 ) ( 32+x2 ) 


32-l 


Taking  the  imaginary  parts  of  either  side,  we  get 


{x  arctan  ( -x/3 ) } /{ ( 1+x2 ) ( 32+x2 ) }dx  = 


2Wc2iv2U^v  -  _ H _  log  ----- 

32-l 


-  00 


23 


- 


' 


' 
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Thus , 


1(6)  =  2 


IT 


1 


6-1 


log  fr dB 


=  2tt 


—  (1-— )_1log(l+i  )d6  -TT  log  2  log 
32  32  6  6+1 


=  2tt(-  +  — 


23' 


66' 


^  +  -1 


36 


203' 


23 

903( 


)  -  TT  log  2  log 


6-1 

6+1 


1  -1 

by  using  the  power  series  expansions  for  both  (l- —  )  and  log  (1+1/6) 

62 

and  the  fact  that  the  right-hand  side  vanishes  as  6"*”.  Specially  if 
6=3,  by  direct  computation 

1(3)  =  1,17^9 
Substituting  in  (*+.2.8), 


Lim  e 
m-*» 


w,w*  (Cauchy)  =  .*+58. 


It  is  easily  seen  that  for  m=2 ,  e^  ^  (Cauchy)  =  ,75.  Thus 
the  efficiency  of  the  joint  ranking  procedure  relative  to  the  separate 
ranking  one  decreases  from  .75  to  .*+58  as  m  increases  from  2  to  infinity  - 
presumably  staying  below  unity  in  the  entire  range. 

Remark  Thus  joint  ranking  procedure,  with  alignment  on  the  mean 
is  very  undesirable  for  Cauchy  distribution.  The  deficiency  perhaps  lies 
in  the  choice  of  arithmetic  mean  as  our  "norm"  of  alignment,  rather  than 
in  the  joint  ranking  procedure. 

The  median  is  a  more  stable  measure  of  central  tendency  for 
the  Cauchy  distribution  than  the  arithmetic  mean  and  may  be  preferable. 

But  the  investigations  of  these  questions  seem  rather  complex  and  will 
not  be  taken  up  here. 


II .  Partially  Balanced  Incomplete  Block  Design 


In  this  case,  the  corresponding  separate  ranking  statistic 


, 


-  4l  - 


is  given  by 


(4,2.9)  W*  =  {I2(k-l)/nm(ni-l) }  l  (R*  -  ^-(m+l)}2 

n  j=i  J 

where  R.  =  sum  of  the  ranks  under  the  j  treatment.  Its  large  sample 
J 

behaviour  is  given  by  the  following  theorem. 

Theorem  H.  Assuming  for  each  fixed  n,  the  truth  of  the  hypothesis 

Kn »  if  (i)  F(x)  possesses  a  continuous  derivative  a,e.  and  (ii)  there 

exists  a  function  b(x)  such  that  [{F(x+0)  -  F(x)}/0]  <_  b(x)  where 
r°° 

b(x)dF(x)  <  °°,  then  the  statistic  W*  given  by  (4,2.9),  converges, 

i  ^ 


as  n-*»,  to  x,2  -.(A*  «),  where 
k-1  w*  5 


(4,2.10)  =  (l2m(m-l) (k-l)/(m+l)k2}  [ 


w 


k 


{  F  ’  ( x )  }  2  dx  ]  2  l  {0.-  /k-r-v 
A  j  m(m-l) 


j=l 


,  1  p . . , 0  .  ,  }2 

L  JJ  J 

j'+j=1 

Proof;  There  are  only  a  finite  number  (m)  ways  of  choosing  m 
treatments  from  k  and  the  vectors  of  observed  ranks  containing;  the  same 
sets  of  treatments  are  independently  identically  distributed.  The 
conclusion  of  the  theorem  now  follows  from  the  central  limit  theorem  for 
random  vectors. 

From  (4. 2, 10)  and  (3,3.3),  taking  the  ratio  of  the  non¬ 
centrality  parameters,  the  efficiency  expression  is  seen  to  agree  with 
that  given  by  (4.2.3),  So  the  comments  of  the  earlier  section  apply 
verbatim  to  this  case. 

The  asymptotic  efficiency  relative  to  the  F-test,  however, 
depends  upon  the  0's  and  thus  does  not  convey  much  meaning. 


■ 


■ 
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HI”  "Balanced"  Incomplete  Blocks  with  p  Observations  per  Cello  (For  the 
meaning  of  "balanced"  see  section  3.4) 

In  this  case  the  separate  ranking  statistic  can  be  written 

in  the  form 


(4.2.11) 


k 

W*  =  {  12/pkp2(mp+l)}  l  {R? 
^  i  J 


rp(mp+l )  ,2 
2  J 


If  the  conditions  of  Theorem  G  are  satisfied,  the  statistic  W*  tends, 

’  n 

as  n-*>°,  to  a  x2^  -^(A"^)  variable,  where 


(4.2.12) 


=  (  12p2m(m-l)/(k-l) (mp+l) } 


(f'(x)  }2dx  |  l  (e.-e) 
Jj=i  J 


2 


This  result  can  be  proved  in  exactly  the  same  way  as  Theorem  G. 

The  classical  F-statistic  is  distributed  in  the  limit,  as 
n-*»,  as  x2^  1(A"f)»  where 

k 

(U.2.13)  A"  =[p(m-l )  T  ( 8 ,-e ) 2 ] /[ (k-l)a2 ] 

F  J=1  J 

if  the  variance  a2  for  the  parent  distribution  exists. 

From  (4.2.12),  (4.2.13)  and  (3=4.3),  taking  the  ratio  of  the 
non-centrality  parameters,  we  get 


(4.2.14) 


e  * 
w,w* 


[  (m-l)(mp+l)  (kp-l) /4mp(mp-l)  (k-l)  (l-3l)  ] 


.[ 


CO 

{G T (x) }2dx/ 

CO 


CO 

(F ’ (x) }2dx]2 

J-  00 


(4.2.15) 


[ 3a2 (kp-l) (m-l) / (mp-1 ) ( k-l) ( 1-3A ) [ 


r°° 

(G  * (x) }2dx] 2 

oo 


These  expressions  are  again  of  the  same  type  as  (4.2.3)  and  (4.2.4). 


It  is  interesting  to  note  that  if  m  and  k  are  kept  fixed 
and  p  allowed  to  tend  to  infinity,  then  using  results  of  Theorem  4.1 
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(4.2.16) 


Lim 

p-*» 


w  ,w 


* 


k(m-l)  ,  , 

5TkHT- 


if  the  parent  distribution  has  a  finite  expectation.  We  also  note  that 

the  efficiency  expression  e  ^  of  (4.2.15),  differs  from  that  of  (4.2.4) 

w ,  r 

by  a  constant  factor  { (kp-l) (m-l)/(mp-l) (k-l) }  which  is  always  less 

than  unity  and  tends  to  (k(m-l) /m(k-l) }  as  p  tends  to  infinity. 

Remark .  The  above  ’’balanced"  incomplete  block  design  with  p 

observations  per  cell  can  be  treated  as  a  partially  balanced  incomplete 

block  with  two  associate  classes  where  p..,  =  r  or  p .  We  have  the 

JJ 

additional  information  that  the  treatment  effects  for  plots  occurring 

in  the  same  cell  are  equal.  However  it  follows  from  (4.2.16)  that  this 

additional  information  regarding  0.'s  contributes  more  towards  the 

J 

efficiency  of  the  separate  ranking  than  the  joint  ranking  procedure. 
Accordingly,  comparing  the  results  of  this  section  with  those  of  the 
complete  case  (Mehra-Sarangi  [l4]),  it  seems  that  joint  ranking  pro¬ 
cedure  (taking  account  of  interblock  comparisons)  would  be  preferable 
to  separate  ranking  procedure  only  if  the  number  of  observations  in 
the  different  cells  are  of  the  same  order. 

We,  however,  wish  to  remark  that  the  limiting  value  of 
e  j,  as  p  tends  to  infinity,  will  be  rather  close  to  unity  if  the 
value  of  m  is  close  to  k. 


4 . 3  Joint  Ranking  Within  Similar  Blocks  in  a  Special  Case 

In  this  section,  two  blocks  are  called  similar  if  they 
contain  the  same  set  of  treatments.  Here,  we  consider  a  method  of 
ranking  in  which  all  aligned  observations  within  each  set  of  similar 
blocks  are  ranked  jointly  together.  Different  similar  sets  are  ranked 
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separately.  A  test  statistic  can  be  proposed  by  finding  the  limiting 

distribution  of  (Ri ,^2 s » * . ,R, )  where  R.  refers  to  the  sum  of  the  ranks 

^  J 

t  h. 

under  the  j  treatment »  We  specially  consider  a  balanced  incomplete 

block  design  with  n=Nd  blocks  containing  d  sets,  each  of  N  similar 

blocks.  Suppose  d  is  kept  fixed  and  n  tends  to  infinity  with  N.  In  the 

notations  of  sec.  3.2,  we  prove  the  following  theorem. 

Theorem  4.2.  Suppose  the  hypothesis  k  Ue  true  and  the 

-  n 

conditions  ( i ) ,  ( i i )  and  ( i i i  )  of  Theorem  3 . 1  be  satisfied .  (See  p . 18 ) . 
Then  the  statistic  V  ,  given  by 


(4.3.1) 


V  =  {3d2(k-l)/m2(m-l)(l-3A)n3}  £  {R .  -  —  (Nm+l)}2 

n  j=l  J 


converges  in  distribution,  as  n-*=°,  to  x2^  -^(A^)  w^ere 


(4.3.2) 


A  =  [ 3(m-l) / (l-3A ) ( k-l) ]  [ 
w 


{ G ' (u) }2du]2  l(B.-B)2 


>1  J 


is  the  same  as  given  by  (3.2.3). 

Proof ;  Since  the  ranking  involves  d  sets,  each  containing  m  out 

of  the  k  treatments  to  be  compared;  each  R.  can  be  regarded  as  the  sum 

J 

of  (md/k)  independent  r.v  '  s.  Using  results  of  Lemma  3=1  (MeHra- 
Sarangi  [l4]),  it  is  verified  by  direct  computation  that  if 


(4.3.3) 

T  .  = 
nJ 

(4.3.4) 

= 

then  T  =  (T  .. 
jp  nl 

,Tn2,e 

mean  vector  y 

j  2 


{ G f (u) }2du]  (0.-0) 
.  J 


'nk' 


vector  y  =  ( y .  ,y „  , .  .  .  ,yn  )  and  variance-covariance  matrix  A  = 
i  z  k 


1 


j  y 


The  rest  of  the  theorem  follows  from  here,  using 


*  ' 
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Lemma  3.2  (Mehra-Sarangi  [lU]),  just  as  Theorem  3.1  followed  from 
Lemmas  3.1  and  3.2. 

It  is  seen  from  (1.3.2)  and  (3.2.3)  that  in  our  special  case, 
joint  ranking  within  similar  blocks  is  asymptotically  as  efficient  as  the 
overall  joint  ranking  procedure.  This  result  is  striking  and  suggests 
investigation  of  the  fact:  does  the  asymptotic  relative  efficiency  of 
the  joint  ranking  procedure  within  similar  blocks  compared  to  the 
overall  joint  ranking  procedure  tend,  as  the  sequence  of  alternatives 
approach  the  null  hypothesis,  to  unity  from  above  or  from  below?  It 
may  also  be  of  interest  to  try  this  method  of  alignment  on  other 
experimental  designs.  But  these  questions  will  not  be  considered  here 
any  further. 
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CHAPTER  V 

CONCLUDING  REMARKS 

5 . 1  A  Comparison  With  Competing  Test  Procedures 

It  is  a  problem  of  vital  importance  to  compare  the  present 
test  procedure  with  other  competing  test  procedures.,  This  question  was 
partially  answered  in  the  last  chapter  by  comparison  with  separate  rank¬ 
ing  procedure  and  F-test,  In  the  present  work  the  efficiency  expressions 
e  r(F)  and  e  *(F)  were  evaluated  only  for  the  case  of  a  normal  parent 
distribution,  The  problem  of  evaluating  them  for  other  parent  distribu¬ 
tions  e.g.  the  uniform,  logistic,  (double)  exponential,  Cauchy  and  others 
needs  additional  investigation .  These  values  are  available  for  m=2  and 
when  m  tends  to  infinity  (Mehra-Sarangi  [ll+]).  The  value  of  e  -  can  be 

W ,  r 

arbitrarily  high  and  even  infinite  as  in  the  case  of  a  parent  Cauchy 
distribution o  It  may  be  possible  to  find  a  lower  bound  for  e  comparable 

W ,  r 

to  the  one  obtained  by  Hodges-Lehmann  in  [7],  for  the  efficiency  of  the 
Wilcoxon  (Kruskal -Wallis )  test  relative  to  the  t-test  (F-test),  The 
Lemma  5»1  below  provides  us  with  a  lower  bound  for  e  p(F)  for  the  case 
when  F  is  a  symmetric  distribution. 

Lemma  5  - 1  If  the  parent  distribution  F(x)  is  symmetrical ,  then 

A  >  1/6  and  e  r(F)  >  0,432  m/m-1. 

—  w,r  — 

Proof,  If  the  parent  distribution  is  symmetrical,  then  both  G(x) 
and  H(x,y)  are  symmetrical  about. the  origin  and  the  line  y  =  x,  thus 
giving 

(5,1=1)  G(x)  +  G(-x)  =  1 


(5.1,2) 


H( -x ,  -y)  =  1  -  G(x)  -  G(y)  +  H(x,y) 
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Hence , 


A  = 


'-00  j  —  00 


G(x)G(y)  dH(x,y) 


o  Jo 


{ G ( x ) G ( y ) +G ( -x ) G ( -y ) }  dH(x,y)  +  |  j  {G(x)G( -y )+G( -x)G(y ) }dH( -x,y ) 


0  J  0 


r  oo 


(l-G(x)-G(y)+2G(x)G(y) }dH(x,y)  + 


0  0 


{G(x)+G(y)-2G(x)G(y)}dH(-x,y) 


0  0 


r°o 


{1  -  2G(x)  -  2G(y)  +  l+G(x)G(y )  }dH(x,y ) 


0  0 


+ 

. 

{CO 

0  J 

•00 

(G(x)  -  G2(x)}dH(x,y)  + 

CO 

{CO 

'■  J 

oo 

+ 

. 

■  00 

0  ^ 

•00 

{G(x)  -  G(y) }2dH(x,y)  + 

o  J 

•  CO 

0  J 

(G(y)  -  G2(y)}dH(x,y) 


{G(x)  -  G(y)}2dH(-x,y) 


1 


OO  { 00 


{CO  {CO 


{2G(x)-l}  {2G(y)-l}  dH(x,y)  +  (G(x)  -  G(y) }2dH(x,y ) 


0  0 


;  0  J  0 


/•OO  {CO 

+  (G(x)  -  G(y) }2dH(-x,y)  >  i 

J  „  J  „  ~  b 


0  0 


since  each  term  in  the  above  expression  is  non  negative.  In  the  above 
simplifications  the  first  integral  sign  contains  the  limits  for  x  and 
the  second  contains  the  limits  for  y. 

Again  , 

r  i  r  ,oo  -2 


3a2/ (1-3A ) 


(G ' ( x ) } 2  dx 

j—  oo 


=  12a2(m-l) { 


^G'  (x)j2dx}2/m  •  {m/(m-l)  }{l/4(  1-3 A )  } 


>  0,1+32  m/(m-l)  . 
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since  the  variance  of  the  r.v,  with  c.d.f.  G(x)  is  (m-l)a2/m  and  thus 
by  the  results  proved  in  Hodges-Lehmann  [T],  the  quantity  within  the 
first  square  brackets  is  not  less  than  0.86U;  and  A  >_  1/6.  The  result 
trivially  holds  if  the  parent  distribution  does  not  have  a  finite 
variance.  This  proves  Lemma  5°lo  It  may  be  possible  to  improve  upon 
the  lower  bound  (as  a  function  of  m)  of  A  by  finding  lower  bounds  of 
the  non-negative  terms  that  were  ignored  in  proving  A  >_  1/6. 

For  the  case  of  normal  parent  distribution,  the  conditional 
test  considered  here  is  slightly  less  efficient  than  the  test  proposed 
by  Lehmann  in  [12]  but  is  more  efficient  than  his  test  proposed  in  [ll]. 
But  it  is  not  at  all  clear  how  the  present  test  stands  in  relation  to 
these  for  other  forms  of  the  parent  distribution.  However  this  condi¬ 
tional  test  has  the  important  advantage  of  providing  exact  levels  of 
significance  for  all  sample  sizes,  whereas  the  tests  proposed  by 
Lehmann  in  [ll]  and  [ 12 ]  are  only  asymptotically  distribut ion- free . 


■ 

a 
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